Topics in non-perturbative QCD 



Adriano Di Giacomo 

Dipartimento di Fisica E. Fermi, Universita di Pisa, and I. N.F.N. Via Buonarroti 2, 

56100 Pisa 

Received XXX 

1 QCD. Perturbative and non perturbative. 

1.1 Basic field theory. 

A quantum system is denned by the canonical variables q(t), p(t), and by the 
Hamiltonian H(p, q). 

q(t), p(t) is a short notation for qi(t), Pi(t), the index i running over the degrees 
of freedom. In a field theory 

qi(t) = (p a (x,t) (1) 

or i = (a, x). There are infinitely many degrees of freedom. 

Solving the system means to construct a Hilbert space on which q, p act as 
operators obeying canonical commutation relations at equal time and the equations 
of motion. A ground state must exist. 

The usual approach in quantum field theory is to split the Hamiltonian H in a 
free Hamiltonian Hq and an interaction Hi 

H = H +H! (2) 

and to assume as Hilbert space the Hilbert space of the free system Ho, and as 
a ground state the corresponding Fock vacuum. Scattering amplitudes between 
fundamental excitations are then computed by perturbation theory. This approach 
is known as perturbative quantization, and is what is presented in textbooks. 

It is imediately apparent that QCD will present problems in this procedure: free 
quarks and gluons are not observed in nature. Perturbing in their Fock space does 
not seem a sensible procedure. 

A more fundamental approach is needed. 

1.2 Feynman-path formulation. |Q] 

The basic quantity in Feynman's formulation is the action 

S[q] = J dtC[q,q] (3) 

where C = pq — H is the Lagrangian. 
S[q] is a functional of q(t). 
In terms of S the equations of motion are 



SS =0 (4) 



5qi(t) 
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The symbol in eq.Q) means functional derivative. 

The solution of the system is equivalent to compute the partition function Z 



Z = 



JJd?i(*) 



exp{i5[q]} 



(5) 



or better the generating functional 



m-J 



i,t 



exp{iS[q]+i / dtJi(t)q z (t)} 



(6) 



with Ji(t) arbitrary functions of time. 

The integrals in (^J), (^) are on a continuous infinity of variables (functional 
integrals). 

Z[J], if it exists, provides the solution of the system, in the sense that all corre- 
lation functions of the fields are known in terms of Z[J] 



1 5 n Z[J] 

Z K ' SJ il (h)...6J in (t n ) 



= (0\T(q il (t 1 )...q in (t n ))\0) 



(7) 



j=o 



The knowledge of the correlation functions is indeed equivalent to the knowledge 
of the Hilbert space, field operators and ground state (reconstruction theorem). 
The proof of (@) goes as follows. 

The correlations functions obey the following set of coupled differential equa- 
tions, with boundary conditions fixed by the T prescription 

2? t (0|r( ? (i) ftl (t 1 )...g in (t B ))|0) (8) 

n 

= ^6(t- t k )Su„ (0\T (q(t)q n (h) . . . q^ (t fc _i)« ifc+1 (t k+ i)q ln (*„)) |0) 



fc=l 



Here 



V t q{t) = 



(9) 



indicate the equations of motion, and putting T> t out of the T product means that 
time derivatives contained in V t must act on the T product. By use of eq.(g) the 
only terms which survive are those obtained by differentiating the 9 functions in 
the T product, giving the S functions in the right hand side, and canonical equal 
time commutators. 

The left hand side of eq. ([?]) has the T prescription built in, and obeys the same 
differential equations, which proves the equality. 
Indeed 



- / [JJ dq}V tqi (t) qil (h) ...q in (t n )exp[-iS(q)] 
1 f 5 S 

Z J [n^5 g Tmfc(*i)---5in(*™) ex P[- i5 '(9)] = 



A 
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J IN {s^ exp[ ' lS{q)] } 9n(<l) • • • qiM = (10) 

/ ill ^ X] ~~ i fc)^ fc exp[-iS , (g)]g ll (ti) . . . g ifc _ I (tfc_i)t?i fc+I (*fc+i)&„ (*n) 
^ fe=i 

The first equality follows from eq.(Q), the third by partial integration. The finite 
term drops if exp[— iS(q)] vanishes at large q's. In the euclidean this amounts to 
have S(q) — > oo as q — > ±00 (stability). The argument is not true for a ip 3 theory 
or if in a ip 4 theory the coefficient of ip 4 has the wrong sign. 

A more direct proof is the original construction of FeynmanjjJ. q is a complete 
set of operators. The amplitude 

t>(q'\q)t = (q'\cM-iH(t' -t)]q) (11) 

contains all the physics of the system. 

We now split the time t' — t into n + 1 intervals of length S, with the idea of 
sending n — > 00 

J = ^ 

71+1 n— >oo 

Then, using completeness 

t '(?'k) t = | dgi...dg«(9 , |e- <H5 |?n><9n|e-" ra |fo-i)...<9i|e- <Ha |«) (12) 
Consider, for the sake of simplicity, p independent forces 



Then 



(q k+ l\e- iHS \q k ) % e^ 5 (q k+1 \e^ 2s \q k ) 



exp < iS 



(gfe+i - gfc) 
2d 2 



2 



-V(g) 



= exp{i<5£[%]} (13) 



and 



~ n 

v {q'\q) t = lim / T[dq k e iS ^ (14) 

J 1 

Eq.(pT[) is the definition of a functional integral as a limit of ordinary integrals. 

The above construction also provides a continuation to Euclidean time. 

Consider the amplitude irW\Q)—iT a s T — ► 00, in the presence of an external 
source J(r) which is non zero in the interval — T/2, T/2. Let \E n ) be a complete 
set of states of given energies 



T lim iT {q'\qU T = ]T (gV^KW' I j [Q d?] e - SW -/ J(t) * Wlft |^)<£?„|e- flT / a | g ) 
(?'|0)(0|g)<0| f [Jjdg] e - s W-/ J W'W««|o) (15) 
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Only the ground state survives the limit. 

1.3 Feynman integral and perturbation theory. 

To describe perturbation theory in the language of Feynman integral, the La- 
grangian is split into a term bilinear in the fields Co plus the rest 

£ = £ + £ I 



or 



where, in generic notation 



S[q] =S Q + S! 



So = :J»Jtf> 



(16) 



(17) 



2 ot/3 

and Si is of higher degree in the fields. The matrix in eq.(fl7|)is denoted M _1 for 



convenience. 
Then 



exp[-5] = exp(-S ) \ ]T 



(~l)" 5 n 



(18) 



A power exansion is performed in Si. The partition function becomes 



Z = J [Y[dq] exp(-S ){£^S? 



(19) 



and 



•*»(*».)) I°> = | J [Y[dq\ e- s ° qil (ti 



(20) 

If 5/ is a polynomial in the fields the generic integral in ( p0[ ) is gaussian, i.e. has 
the form 

[l[dq i e-^ M ^q ai ...q an = ^ J] M« iQj (21) 



The sum is extended to all possible set of pairs of the indices, M is the propagator, 
and ( pl| ) is nothing but Wick's theorem. 

A technical point arises in gauge theories where the matrix M _1 is not invertible, 
and the naive propagator does not exist due to gauge invariance: the way around 
is Faddeev-Popov procedure 0], which brings back to eq.(|2l]), with a well defined 
M, and the addition of ghost fields. 

The procedure explained above, to expand the weight function in Z in powers of 
Si is far from harmless. This is physically expected, because perturbing around the 
Fock vacuum of quarks and gluons and computing scattering amplitudes thereof, 
is not a good approximation in a theory which confines them. 



A 
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This desease manifests itself in the fact that renormalized expansion of a generic 
observable O in powers of the coupling constant a s = g 2 /47r 

Oj2<Or (22) 
is not convergent, not even as an asymptotic series 

1.4 Feynman path and lattice QCD. 

The lattice formulation |4| is a discrete approximant of QCD Feynman integral, 
in the sense ofeq.'s (H), @. The approximant is expected to be good if the lattice 
is large compared to the correlation length A and A ^> a, the lattice spacing. The 
building block of the theory is the parallel transport between two neighbouring 
sites of the lattice 

U^{n) = exp [iagA^n)] (23) 
where ft is the site, ji = 1 . . . 4 the direction to the next site, and 

a 

with T a the generators of the group in the fundamental representation. The ac- 
tion is related to the parallel transport around the elementary square surface Il^u 
(plaquette) 

2 N 

S = (3 Tr{l-IV} P=-£ (24) 



r 



In the limit a — > 



S ^n-T E G^G^ + Oia 6 ) (25) 



a^Q 4 



i.e. S differs from the continuum action by "irrelevant" terms. 

It can be shown that the lattice action defines a theory with an hermitian Hamil- 
tonian, or that, in the language of statistical mechanics, a transfer matrix exists. 

The theory has a mass gap: the string tension a is observed in numerical simu- 
lations. 

The theory has a fixed point as (3 — > oo, where the correlation length diverges. 
In formulae the lattice spacing is related to the physical scale A as 

a = jcxp{-b {3) (26) 

where — 60 < is the first coefficient of the pcrturbativc f3 function. 

Eq. (|6|) means that as the critical point (3 = 00 is approached the lattice spacing 
becomes exponentially small in physical units, and the coarse structure of the lattice 
disappears (scaling region). 

In summary if (3 is sufficiently large, so that the physical correlation length A is 
such that A > a and at the same time A ^ L, the size of the lattice, the lattice 
should provide a good approximant to continuum Feynman integral. 
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The partition function 

Z = J [n^(«)]ex P [-5] (27) 

is finite, since the group is compact. 

The formulation is gauge invariant and needs no Faddeev Popov. 

Of course it would be perfect if anybody could compute Z analitically, and then 
perform the continuum limit. The problem is open for competition. 

A practical approach Q is to compute Z numerically, by producing via Mon- 
tecarlo techniques a "significant sample" of configurations on which expectation 
values can be computed. As we shall see below if properly asked lattice can give 
important answers. 



1.5 Non convergence of perturbative expansion. Non perturbative formula- 
tion of QCD. 

As first noticed by Dyson |6| QED cannot be analytic in a at a — 0. A gas of N 
electrons in a given volume has energy 

E oc cN + daN 2 

with c, d positive coefficients. The second term stems from Coulomb repulsion. If 
a —>■ —a, the system becomes unstable, and hence no neighbour of the point a = 
exists in which physics is analytic in a. 
Still QED works as an effective theory. 

As we will see below a similar phenomenon occurs in QCD. Perturbation theory 
is ill-defined. Still it works at short distances, for reasons which are not really 
understood. In addition QCD, as argued in sect. 4, is not an effective theory, but 
exists in the sense of constructive field theory. 

Lattice is the most rigorous approach beyond perturbation theory. There exist 
other attempts in the same direction, which we will briefly comment upon 

1. A phenomenological approach is the SVZ sum rules (tJ: non perturbative ef- 
fects are parametrized by condensates, which are defined via Wilson Operator 
Product Expansion (see sect. 2.1 below). 

2. An interesting approach is known as iV c — > ooQ. The idea is that the structure 
of QCD is almost independent on the number of colours N c : the theory at 
N c = oo differs from that at finite N c by small corrections C(l/iV c ) which 
do not alter the main physical features. More precisely the limit is defined as 
N c — > oo at g N c ~ A fixed. Only planar graphs survive at N c — oo. Quarks 
loops are non leading in this expansion. Indeed a quark loop with k gluon 
vertices is 0(N f /Nc /2 ). 

This idea is supported by lattice simulations: apart from a rescaling due to 
the difference in the beta function quenched simulations (no quarks loops) 
differ from full QCD simulations tipically by 10%. 



G 



A 
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The N c — > oo idea has lead to the solution of the U(l) problem ||, [n|. Al- 
though it is difficult to perform calculations in the limt, a number of very im- 
portant qualitative results can be obtained in this approach, which go deeply 
into the structure of the theory. 

A number of non perturbative models exist of QCD which assume that some 
amplitudes are dominant and provide a good approximation to the full theory. 

Among them we quote the instanton gas (or liquid) modeljll]] , and the stochas- 
tic vacuum model |]l2|| . 

3. Instanton gas (liquid) model. Perturbation theory (sect. 1.3) can be viewed 
as a saddle point approximation to the Feynman integral, approximating it 
with a gaussian expansion around the trivial saddle point with zero fields. 
For memory the saddle point approximation to an integral is 



dxg(x)e iJ W = / dSg(x + S)e^ f ^ + ^ 



E 



(iA/)' 



(28) 



with 



Af = f(x + 8)-f(x) 



2 



and x a saddle point f'(x) = 0. It reduces to a gaussian integral if g and Af 
are expanded in powers of 8. 

In principle all saddle points should be included, and their contribution added 
to get a better approximation. 



For the Feynman integral eq.(28) reads 

Z = J [dq]e- S ^ : 
where 



jexp 



2 SqiSqj 



E 



i-iy 



(AS) n (29) 



ss 

~5q 



= 



S = S(q) 



1 r? 2 9 

AS = S[q + 8} - S - --^-^8,8 



q=q 



2 SqiSqj 



) 3 



Extending the approximation to all saddle points means to go beyond per- 
turbation theory: all the solutions q of the classical equations of motion 

^ =0 have to be included, which have finite action S. If S — oo 

the factor e~ s in ( p9| ) will kill the contribution. 

Classical solutions with finite euclidean action are called instantons. they 
have non trivial topology. Indeed finite action means finite J F 2 d i x. To have 
convergence 



F 



O 



1 
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As a consequence A p is a pure gauge on the sphere S3 at infinity, resulting in 
a mapping of S3 on the gauge group. The jacobian of this mapping is given 
by 

dU _ 

where dU is the volume element of the group, normalized to J dU = 1, da^ 
the element of the surface. The chern number 



S3 



k^da^ = n 



is an integer (topological charge) and counts how many times the group is 
swept in the mapping. By Stokes theorem 



k ft da^= (d^k li )dV = n 
s 3 J 

Q(x) — d^k p is the density of topological charge. In QCD 

Q(x) = ^ G %G% (30) 

where = \e[ivpaG a prj is the dual of G a pv . 

Topology ensures stability of the classical solutions. 

Explicit solutions (B.P.S.T. instantons) were first found in ref.Q. 

A line of research attempting to classify all possible instantons developed 
immediately after the appearance of ref. [^3|. A model was strongly pushed 
assuming that the Feynman integral is dominated by quasi solutions, with 
instantons at distances large enough that they can be considered as indepen- 
dent (instanton gas)|0. The attempt was frustrated by infrared problems, 
and by the failure to describe confinement. The model was subsequently re- 
placed by allowing overlapping between instantons (instanton liquid), which 
is useful to describe chiral properties of the theory[pr|. 

Stochastic vacuum. The idea is to compute physical quantities in terms of 
gauge invariant field strength correlators. The simplest of them is 

2V, PCT = (G fW (x)Sf(x,0)G pr7 (0)) (31) 

whith Sq^ the parallel transport in the adjoint representation from to x 
along the path C . A cluster expansion is then performed, and higher irre- 
ducible correlators are neglected. 

Correlators extracted from lattice simulations p5| are used as an input. 

A good phenomenological model to describe high energy scattering between 
hadrons and heavy quarks physics Q. 
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2 Non perturbative effects in QCD. Vacuum condensates. 
2.1 The SVZ sum rules. 

The vacuum condensates were first introduced in rcf. Q to parametrize non per- 
turbative effects. 

Consider the Wilson operator product expansion (OPE) of the product of two 
e.m. currents 

TV»(x)r(0)) ~ Cr(x)-I + C£ v (x)^G^G« a + Cg(x)mW + ... (32) 
\x\->o g w 

The product is written as a sum of local operators, ordered by increasing dimension 
in mass, multiplied by coefficient functions. 

Wilson OPE is a theorem in perturbation theory, and is presumably valid in 
QCD at short distances. 

In Fourier transform, taking the vacuum expectation value (vev) and renormal- 
izing one gets 

IVte) = J d'xe^iOlTU^MO)) |0) = (q^q v - q 2 g^)fl(q 2 ) (33) 

Here n(g 2 ) = n(g 2 ) - 11(0) and, from Eq.@ 

n( (7 2 ) = C 1 (0|/|0) + C G ^ + C v; ^ (34) 

Ci, Cg, C^, are dimensionless coefficients, which are to be computed in perturba- 
tion theory and (0|/|0) = I, 

G 2 ^ {0l M :G a G a \ 0) (35) 

9 

has dimension 4, is the vev of the dilatation anomaly and is known as gluon con- 
densate. 

Gj, = (0|m : : |0) (36) 

has also dimension 4, is related to the spontaneous breaking of chiral symmetry, 
and is known as quark condensate. 

G2 and G^ are not defined in perturbation theory, where are usually put equal 
to zero, so that only the first term survives in eq.(|34]). They parametrize non per- 
turbative effects. 

Eq.(|3) is the basis of SVZ sum rules. 

The l.h.s. is parametrized by a dispersion representation 

urt . v r * (37) 

J 4ml s s — Q + ie 

where 

fl( s ) = ^e-^his) (38) 
^e+e-^^+p- (S) 
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In perturbation theory R(s) is a constant, modulo log's, a oc 1/s, there is no scale, 
and in the r.h.s. of eq.(|34|) only C\ is present. In nature scale invariance is broken 
by resonances: in the r.h.s. this reflects in the presence of condensates. 
An appropriate weighting of both sides of the equation (34) 

2— ^ Ci + Cg-t + C^-j (39) 

4m 2 s s-q 2 +ie q 4 vv q 4 

which emphasizes the ranges of q 2 in which both descriptions are valid, relates 
the parameters of the resonances to the condensates. A successfull phenomenology 
emerges, and a determination of G2, G^ [p7| 

G 2 = (.024 ± .011) GeV 4 (ipij)) = -.13GeV 3 (40) 

In fact the calculation of C\ is affected by the presence of renormalons, an ar- 
bitrariness which mimics the condensates, which, as a consequence, are not well 
defined. 

In spite of that sum rules work and give a consistent determination of the con- 
densates, independent of the choice of the currents considered. 

2.2 Determination of G2 on the lattice. 

Consider the gauge invariant field strength correlators 

V^ pa {x) = (0|Tr \ [ G^{x)Sc(x,0)G pa (0)Sl(x,0)\ |0) (41) 

where 

a a 

with T a the generators in the fundamental representation and Sc a Schwinger 
parallel transport from to a; along the path C 



S c (x,0)=Pexp(i J 



A p {x)dx» ) (42) 

c 



In general the definition depends on the choice on the path C . This will not be 
relevant for the following. For C we shall assume a straight line. T> liv , pcr (x) can be 
viewed as a split point regulator a la Schwinger of G2. The OPE at small x gives 



(x) ~ — (I) + c> (G 2 ) + • ■ • (43) 



with c, d computable coefficients in perturbation theory. (G2) is defined if c, d 
can be unambigously computed. We will see that this is not the case. A generic 
parametrization dictated by symmetry is |]l2|| 

ZW«rO»0 = (9^9 P a - 9^ P )[D(x 2 ) + D^x 2 )] (44) 

dDi 
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Choose x^ along axis. Then|Tf 

P||(x 2 ) = -V 0l V 0l =V + V 1 +x z ^- (4:>) 
1 3 

T> ± (x 2 ) = - ^ 23„ D, 3 = 23 + V x 



On the lattice G MJ/ is a plaquette, the parallel transport is provided by the links 
and the lattice version of 23^, pcr of the correlator is given in pictures by 



V 



l 



The term subtracted insures that the quantum number exchanged is in the ad- 
joint representation. 

At large [3 = 2N c /g 2 (see sect. 4) 



K^po =a 4 2W+0(a 6 ) (46) 
. U{3) 

They can be parametrized as 



....,1 T)L !",,,■ „,„„„.!„„! CT 



23m and 23^ for quenched SU(3) are shown in fig.l 



23f = ^e- x ' x ' + cG 2 e~ x ' x 
a; 4 

Fig. 2 shows the behaviour in full QCD, 4 quark species. 
The corresponding results for G 2 and A are 

SU{2) quenched[18] G 2 = (0.33 ± 0.01)GeV 4 A = .16 ± .02 fm (47) 
517(3) quenchedH G 2 = (0.15 ± 0.03)GeV 4 A = .22 ± .02 fm (48) 
SU(3) full QCD||| G 2 = (0.022 ± 0.005)GeV 4 A = .34fm (49) 

The value obtained from SVZ sum rules is|i~7| 

G 2 = (0.024 ±0.011)GeV 4 

in excellent agreement with eq.(^|). 

However, like in sum rules the determination of the coefficients c in eq.(^3|) is 
again ambigous by terms which mimic the second term in (G 2 ), due to renormalons. 
Nevertheless everything works and compares satisfactorily with sum rules. 

An alternative determination can be obtained by measuring the expectation 
value of the plaquette (density of action). Indeed Q 

1 - IV oc g 2 G^G> 4 

so that 

(0|1 - IV|0) oc G 2 a 4 
Czech. J. Phys. 50 (2000) A 11 
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Here again a term exists which scales like a 4 (eq.(p6|)), but a perturbative contri- 
bution, corresponding in the continuum to a quadratic divergence is present, which 
must be subtracted. This term is again ambigous by contributions which mimic G^. 
Also here, however, the determination is in agreement with (f47f), ([48f). 

2.3 Lattice determination of the chiral condensate {tjjtp} 

The gauge invariant quark correlator 

S{x) = ($(x)S(x, 0)^(0)) (50) 

can also be determined on the lattice pl|. 
In the continuum limit 

S{x) ~ 4 + (51) 

On the lattice 

S L ( X ) = + a 3 W)e- x/A > + 0(a 5 ) (52) 



The result of the fit gives 

\(M 1/3 = I 



(205 ± 21)MeV am = 0.02 J = 0.65(3) 
(160 ± 12)MeV am = 0.01 ^ = 0.57(2) 



A = 1.4m" 1 (53) 

A more careful analysis of the chiral limit is needed. 

In conclusion the condensates can be extracted from first principle on the lattice, 
and agree with the phenomenological determinations, in spite of the fact that G2 is 
undefined by renormalon ambiguities both in the lattice determination and in the 
sum rules. 

The lattice correlators provide the gluonic correlation length. In quenched QCD 
they are smaller by a factor 3 with respect to full QCD. 

The lattice correlators are used as an input in the stochastic vacuum model. 

3 Confinement of colour. 

3.1 Introduction. 

Quarks and gluons have never been observed: asymptotic states in QCD are 
colour singlets. This phenomenon is known as confinement of colour. 

In nature the ratio of the abundance of quarks n q compared to that of nucleons 
n p has an upper limit p2| 

^ < 10 - 27 (54) 
corresponding to Millikan analysis of ~ 1 g of matter, in search of fractional charges. 
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In the absence of confinement the Standard Cosmological Model predicts n q /n p ~ 
10~ 12 . The argument is as follows) 23 . Let <tq = lim v ->oVcr, <j being the cross 



section for the process of hadronization in quark gluon plasma qq — ► hadrons, 
qq — ► q + hadrons. Quarks will decouple in the process of cooling of the universe 
when the rate of burning is of the order of the rate of expansion 



since n 7 = T 3 , this gives 



a n q = I 



i q _ G 1 ^ 2 _ 10- 



n 



1 



a T 2 a m p T 



An estimate for oo is oo < m ^ 2 - A conservative value for T is T ~ 10 GeV and 
this give 



In physical units (eq.(|26| 



^ > io- 21 

or, since n p /n 7 ~ 10~ 9 , n q /n p > 10 -12 . 

A factor 10 -15 between expectation and the empirical upper limit is too large 
to be explained by tuning small parameters. It can only be explained in terms of 
symmetry, like e.g. the smallness of resistivity in a superconductor |p4| . 

The deconfined phase is observed on the lattice, by simulating QCD at high 
temperature. 

The partition function is obtained by computing the Euclidean Feynman integal 
at infinit spatial extension, on a time interval 1/T, T being the temperature. 

This is done on a lattice Nt X Ng with Ns ^S> Nt and, in terms of the lattice 
spacing a 

— = N T a 
T 

X = ^ exp(6o/3) (55) 

Due to asymptotic freedom high temperature corresponds to high /3, i.e. to weak 
coupling ((3 — 2N c /g 2 ), low temperature to strong coupling. 

The deconfming transition is detected by measuring J25| either the string tension 
a or the expectation value of the Polyakov line (L) , which is the trace of the par- 
allel transport from — oo to +oo along time axis, closing trough periodic boundary 
conditions. 

The string tension is the constant appearing in the linear static potential re- 
sponsible for confinement, 

V(r) = or (56) 

has dimension m 2 , is related to the slope of the Regge trajectories S (~ 1 GeV 2 ) 
by the relation 

°=h s (57) 
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(T =/= means confinement. 
The Polyakov line 



(L) = (Tr { Pexp 



N T a 

A (x,t)dt 

-N T a 



(58) 



vanishes if Zn is a symmetry. This happens in the confined phase. 

In the confined phase (L), which can be interpreted as exp(— /Lt g ), fi q being the 
chemical potential of a quark, has zero expectation value since an infinite energy is 
required to add a free quark to the system (confinement). 

In the language of statistical mechanics a parameter which is different from 
zero in the weak coupling regime, like (L), is called an order parameter, while a 
parameter like a which is non zero in the strong coupling regime is called a disorder 
parameter. 

The meaning of these denominations will be clear below, when the concept of 
duality will be introduced. 

A drawback of both a and (L) is that they only work for pure gauge theories 
in the absence of quarks. Indeed if dynamical quarks are present the string tension 
is not defined, since taking qq apart does not produce a potential energy but qq 
pairs. With dynamical quarks Zjq is not a simmetry any more, since quarks belong 
to the fundamental representation, which has non zero triality, and cannot be an 
order parameter. A true order disorder transition is related to chiral symmetry, the 
corresponding disorder parameter is (ipip) > but it is not precisely understood what 
chirality has to do with confinement. 

A genuine order (or disorder) parameter for confinement is needed, which can 
work both with and without quarks. 

Indeed, in the spirit of N c — > oo, as described in sect. 1.5, the underlying dynamics 
must be the same, independent of N c and Nf, for sufficiently small Nf, apart from 
small corrections 0(1/N C ). 

This will be one of our main prejudices in the search of a mechanism for con- 
finement. The other prejudice will be that a suppression by a factor 10 -15 of n q /n p 
can only be explained in terms of symmetry. 

The symmetry of the confined (disorderd) phase can be understood in terms of 
duality ||, |7|. 

Duality is a deep concept in statistical mechanics and in field theory. It applies 
to systems which admit non local excitations with non trivial topology. The idea is 
that, besides the usual description in terms of local fields appropriate to the weak 
coupling regime, a dual description is possible, in which topological excitations 
become local, and local fields become topological excitations. In the dual description 
the original strong coupling regime becomes weak coupling. 

We shall illustrate the concept of duality on a simple system, the 2d Ising model, 
in the next section. 

As for QCD if the dual were known, i.e. if the non trivial topological excitations 
were identified, which interact weakly in the confined phase, an effective Lagrangian 
could be written, and the theory would be essentially solved. 
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This already happens in N = 2 supersymmetric version of the theory |2q|. 

A big step forward has been done in that direction in ordinary QCD, on the 
line of understanding the symmetry of the dual excitations |29|. The results will be 
presented in sect. 3. 3 

3.2 Duality: the 2d Ising model. 

To illustrate the strategy used in QCD and the meaning of duality, we shall 
describe in some detail the 2d Ising model, which is the prototype system showing 
dual behavioural, We will do that by a tecnique allowing to explicitely compute 
the disorder parameter interms of spin variables |]30|] . 

The two dimensional Ising model is defined on a square lattice, by a held variable 



cr in 



±1 



The partition function is 
with 



Z = exp[—S[cr]/T] 
S[<t] = J a(n)a(n + fi) 



(59) 

(60) 
(61) 



J > 0, n is the site and /t are the two independent vectors joining a site with the 
neighbouring sites. 

The model is exactly solvable, and presents two phases separated by a critical 
temperature T c . Defining the magnetization as 



lim 

V->oc V 



for T <T C (cr) ^ (ordered phase) 
for T > T c (a) = (disordered phase) 

The critical temperature is related to the scale J of the model by the equation 

2 J 



Tc huS/2 + 1) 

As T approaches the critical temperature T c (a) tends to zero as 

1 



(a) r ^ T (T c -T) 







(3 



The correlation length £(T) defined by the relation 

\fi — m\ 



(a(n)cr(m)) ~ exp 

\fi— m| — >oo 



behaves as 
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(62) 
(63) 

(64) 

(65) 
15 



Adria.no Di Giacomo 



A field theory in 1 + 1 dimensions is defined at the critical point, where the cor- 
relation length diverges (eq.(|6^)), and the coarse structure of the lattice becomes 
irrelevant. 

The action can be written, up to an irrelevant constant as 

S = ^X> p a(n)] 2 (66) 

rt,p, 

with 

A At (r(n) = a (ft + ft) — <j(ft) (67) 
The equation of motion reads 

A M A M cr(n) = (68) 

and the current 

is identically conserved 

= (70) 

The correspnding constant of motion is 



Q = J dxjo(t,x) = a(x — +oo) — a(x = — oo) (71) 

Q is the number of kinks minus the number of antikinks, and has a topological 
meaning. 

The dual description! 27 of the system is constructed as follows. A dual lattice 
is defined by associating a site of it to each link of the original lattice and a dual 
variable cr*(j) on the dual lattice via its correlation functions 

(a*(iytf)> = ! (72) 

Z is obtained from Z, eq.(|60|) by changing the sign of the links [Ja(ft)a(ft + p,)] 
along an arbitrary path in the dual lattice joing i to j. Then |p7f : 

(i) a* is a dicotomic variable like <r, cr*(i) = ±1- 

(ii) Z is independent of the choice of the path from i to j on the dual lattice. 

(iii) In the thermodynamic limit (V — * oo) 

Z[a,T] = Z[a*,T*] (73) 

with 

sinl 4 = ^ (74) 

or 

T ~ (75) 
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The partition function of the new variables has the same dependence on a* as the 
original partition function had on er, except that high temperature (strong coupling) 
is mapped into low temperature and viceversa. 

(a*) =0 T <T C (a*) ^ T > T c 



or 



= 



(76) 



(a*) is the order parameter in the dual description and is called a disorder param- 
eter. Strong coupling in terms of a becomes weak coupling for a* and viceversa. 

An explicit construction of (a*) can be done which evidences its meaning of 
creation operator of a kinkpC[. As a consequence (a*) ^= is nothing but conden- 
sation of kinks, or spontaneous breaking of the topological symmetry fl70|). In the 
hot phase the vacuum is a superposition of states with different values of Q, the 
number of kinks. 

The disordered phase looks ordered in the dual description. 

The operator which creates a kink at x = x and time to can be writtenfeoj 



,(*) 



exp 



2 — AqOMo) 

X>Xq 



(77) 



Ao(x, to) is nothing but the conjugate momentum to a with the Hamiltonian eq.(|6q) 
and ( |77| ) is nothing but the version for a discrete euclidean field of the well known 
relation 

e ipa \x) = \x + a) (78) 

by which the position is shifted. 

In field theory the analog of x is the field $, the analog of p its conjugate 
momentum, and eq. (f78|) is the tool to add a classical topological excitation to the 
field configuration. Some care is necessary when the field is compact or even 
discrete. 

The expectation value of p can be computed numerically: a more convenient 
technique is to compute the quantity 



(79) 



(//) = exp 



(80) 



The phase transition to the ordered phase produces a sharp drop in (/i), or a 
negative peak on p. 

p vs T is shown in fig. 3 for different sizes of the lattice. At large p can be 
computed in perturbation theory, and the comparison with data is shown in the 
figure. As L — > oo, p — * — oo, or (p) — > 0, as expected. At small (3 p is a constant 
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compatible with zero, and lattice size independent, corresponding to (/i) = 1. In 
the transition region where (/x) - (T c - T) s a finite size scaling analysis can be 
performed to extrapolate to the thermodynamic limit 



sS-l L a 



(jj.) = ( j 9 - J 9)VJ- -1 (81) 
If | < 1, / only depends on L/£ = L(l - T/T c ) u , and 

( M >~(/3-/3 c ) 5 /[£ lA/ (&-/?)] (82) 

or 

JL; ft) (83) 

This relation allows to determine both v and f3 c — f3 from data coming from different 
lattice sizes. 

The result for scaling is shown in fig. 4. The result is compatible with the known 
values 

v = 1 f3 c = 0.44068 

and gives 5 = 0.120(5) to be compared with the exact value S ex = 0.125. 

The construction of the disorder parameter by the rule ([78]) has also been suc- 
cessful! for other systems which are expected to have a dual behaviour: the xy model 
in 3d (liquid i/e4^l|), the compact U{1) gauge theory in 4d|32|], the Heisenberg 
model ||. 

In each case the proper topological excitations condense in the disordered phase. 
We have then used the same method for QCD. 



4 Duality and confinement in QCD 

As explained in sect. 3.1, we expect confinement to be produced by a symmetry 
property of the vacuum, related to the condensation of topological excitations. 
Natural topology of three dimensional configurations comes from the mapping of 
the sphere at infinity, S%, on a group. The corresponding topological quantum 
number is magnetic charge. 

Condensation of magnetic charges in the vacuum is the magnetic analog of con- 
densation of electric charges, which is known as superconductivity, and is named 
"dual superconductivity" |3j . 

The picture fits with the original proposal of ref's.|35|, |36), which ascribed con- 
finement to the formation of dual Abrikosov lines, by the chromoelectric field of a 
qq pair, giving an energy proportional to the distance, or a string tension (eq.(|5^)). 

Monopolcs in a non abelian gauge theory exist as stable particles in the Higgs 
phase of a theory coupled to a scalar field in the adjoint representation |37j [58|. In 
the absence of Higgs field a conserved magnetic charge can be defined anyhow by 
the procedure known as abelian projection. We shall briefly summarize it for the 
simple case of SU(2); the extension to higher groups is routine p9|. 



18 



A 



Czech. J. Phys. 50 (2000) 



Topics in non-perturbative QCD 



A gauge invariant field strength F^ y can be defined|37j as follows 

Fpu = - H(D M $ A DA) (84) 

here is the field strength tensor 

G>„ = d^Av - dvAp + gA„ A l v (85) 
$ is any field belonging to the adjoint representation and 

$ 

$ = (86) 
|$| 

is the unit vector of its direction in colour space. $ is defined everywhere except 
where $ = 0. 

D^8^+ gA^ A (87) 

is the covariant derivative. 

Both terms in (§3) are separately gauge invariant. 

The combination is chosen in such a way that terms in A^A V cancel. It is easy 
to check that 

Ffiv = HdpX, - d v A^) - -$(a M $ A dA) 



The dual of F^ 



f) 



Ff_uj — c^fivpoFpa (89) 



defines a magnetic current 

3u = d„F; v (90) 

which is identically conserved. 

The corresponding magnetic charge is a constant in time. This is true for any 
choice of the local operator <&. There is an infinite number of topological constants 
of the motion in QCD. In the gauge in which $ is constant, e.g. oriented along the 
3 axis, the second term in eq.(pS|) drops and F^ becomes an abelian field 

Fpy = d, t Al - d v A\ (91) 

The corresponding gauge transformation is called abelian projection. If $ were a 
Higgs field, the gauge would be the unitary gauge. 

It can be shown that monopoles can exist at the sites where $ = 0. $ = 
describes world lines of monopoles. 

Of course people like to think of monopoles as stable, observable particles. This 
is not the case in general for the monopoles defined by the abelian projection. 
However the corresponding magnetic symmetry exists and is well defined. One can 
then investigate if the vacuum is invariant under it, or it is a superposition of states 
with different magnetic charges, like the Bogoliubov vacuum of superconductivity. 
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This will in any case give information on the symmetry pattern of the confining 
mechanism. 

As already noticed, there are infinitely many magnetic symmetries, one for each 
abelian projection. There has been a tendency, during the years to think that some 
abelian projection is better than others, and really identifies not only a symmetry, 
but also the real excitations which would be weakly coupled in the dual descrip- 
tion of QCDj3{|. The idea was supported by the so called abelian dominance and 
monopole dominance. The U(l) degrees of freedom emerging from this projection 
and in special those related to magnetic charges seem to account of a large fraction 
of observed quantities (string tension, condensates. . . ). 



An alternative idea|3j§ is that all abelian projections define magnetic charges 
which are physically equivalent. This would mean that the yet unknown excitations 
which become the fundamental fields in the dual description of QCD, have non zero 
magnetic charge in all abelian projections. If true this would be a very important 
step in understanding the dual description of QCD. We have constructed [^9) a dis- 
order parameter for the generic magnetic symmetry, defined by abelian projection, 
by using the techniques illustrated in sect. 3. 2. We have then checked its behaviour 
as a function of the temperature by lattice simulations. We have done that for a 
number of different independent choices of the field $ which defines the abelian 
projection, both in SU(2) and in SU{3) pure gauge theory^] and in full QCD]4C|]. 
The behaviour of p across a deconfining transition is shown in fig. 5. 
The finite size analysis (fig. 6) allows to determine the critical index 8 of the dis- 
order parameter, the critical index v of the correlation length and j3 c . The latter two 
quantites are known by independent method^lj. The results can be summarized 
as follows 

(i) The determinations of (3 C and v are consistent with other methods, showing 
that we are really observing the deconfining transition. 

(ii) Condensation exists for all the abelian projections. 

(iii) The critical index 6 is independent of the abelian projection. 

The definition of the disorder parameter is legitimate also in the presence of quarks. 
There is a peak in p, analogous to that of fig. 5 at the chiral phase transition. Analysis 
is in progress to determine the critical index v and to check that f3 c coincides with 
that of the chiral transition. Our disorder parameter describes confinement also in 
the presence of quarks, and this is consistent with the idea of iV c — > oo. 



Confinement is one of the most fascinating problems in theoretical physics. 
Lattice simulations have provided important hints to understand its mechanism. 
The transition confincmcnt-dcconfincmcnt is for sure an order disorder transi- 



Duality seems at work, the confined phase breaks the magnetic symmetry defined 
by all abelian projections. The dual excitations have magnetic charge in all of them. 




5 Conclusions. 



tion. 
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The dream is to identify the dual excitations, as was done in TV = 4 SUSY QCD 
in ref.j28|. This would give the solution of QCD in terms of a weakly interacting 
effective lagrangian of dual fields. 
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Figure Captions 

Fig.l Field strength correlators in pure gauge SU(3), ref . flij . 
Fig. 2 Field strength correlators in full QCD, ref. Jl9| . 
Fig. 3 2d Ising model, p vs T for different lattice sizes. 
Fig. 4 2d Ising model. Finite size scaling, eq. (|85|) . 

Fig. 5 SU(3) gauge theory, p vs j3 for different abelian projections. The peak signals 
the confinement. 

Fig. 6 SU(3) gauge theory. Finite size scaling of p. 
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